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Abstract. Given a simple polygon P consisting of n vertices, we study 
the problem of designing space-efficient algorithms for computing (i) the 
visibility polygon of a point inside P, (ii) the weak visibility polygon of 
a line segment inside P and (iii) the minimum link path between a pair 
of points inside P. For problem (i) two algorithms arc proposed. The 
first one is an in-place algorithm where the input array may be lost. It 
uses only O(l) extra space apart from the input array. The second one 
assumes that the input is given in a read-only array, and it needs 0(y/n) 
extra space. The time complexity of both the algorithms are 0(n). For 
problem (ii), we have assumed that the input polygon is given in a read- 
only array. Our proposed algorithm runs in 0(n 2 ) time using O(l) extra 
space. For problem (iii) the time and space complexities of our proposed 
algorithm are O(kn) and O(l) respectively; k is the length (number of 
links) in a minimum link path between the given pair of points. 

1 Introduction 

Visibility is one of the foundational areas in computational geometry and it 
has applications in various domains, including robot motion planning, guarding 
art galleries, computer graphics, GIS, sensor network. For an illustrated survey, 
see [4]. Recently, visibility algorithms are being embedded in the hardware of 
digital cameras, sensors, etc, and the constraint on the size of the instrument 
has become important. So, the algorithm designers are now becoming interested 
in developing space-efficient algorithms for various visibility problems. 

So far, in-place algorithms have been studied for a very few problems in computa- 
tional geometry (see [3]). In [1], constant work space algorithms for the following 
visibility related problems are studied: (i) triangulation of a simple polygon, (ii) 
triangulation of a point set, (iii) Euclidean shortest path between a pair of points 
inside a simple polygon, and (iv) Euclidean minimum spanning tree, where the 
input is given in a read-only array. The time complexity of the algorithms for 
the problems (i)- (iii) arc 0(n 2 ), and that for problem (iv) is 0(n 3 ). The open 
question was whether one can compute the visibility of a point inside a simple 
polygon in sub-quadratic time, where the polygon is given in a read-only array 
[1]. Recently, two algorithms for this problem are proposed by Barba et al. [2]. 
The first one is deterministic, and it requires 0{nf) time and 0(1) space, where 
f is the number of reflex vertices of the output visibility polygon. The second 
one is a randomized algorithm and it requires 0(n log r) time and O(logr) space, 
where r is the number of reflex vertices in the input polygon. 



New Results: In this paper we present the following results: 

• An in-place algorithm for computing the visibility polygon of a point inside 
a simple polygon P in O(n) time and 0(1) extra work-space, where n is the 
number of vertices of P. Note that, after the execution of the algorithm, the 
polygon can not be retrieved from the content of the array. 

• An 0(y/n) space algorithm for computing the visibility polygon of a point 
inside the polygon P in O(n) time, where the vertices of P are given in a 
read-only array. 

• An 0(n 2 ) time and 0(1) extra work-space algorithm for computing the weak- 
visibility of the polygon P from a line segment inside the polygon, where the 
vertices of P are given in a read-only array. 

• An 0(kn) time and 0(1) extra work-space algorithm for computing a mini- 
mum link path between a pair of points s and t inside P, where the vertices 
of P are given in a read-only array and k is the size of the output. 

2 Visibility of a point inside simple polygon 

Let the vertices of the polygon P = {pi,P2, ■ ■ ■ ,Pn\ be given in an array P in 
anticlockwise order. Initially, P[i] contains the coordinates of the i-th vertex pi 
of P in the given order. Let tt be a given point inside P. We will consider the 
problem of computing the visibility polygon of tt in P. It is easy to see that the 
visibility polygon of tt can be computed in 0(n 2 ) time and O(l) extra space 
where the vertices of P are given in a read-only array. We will describe two 
algorithms, namely InplaceJVisibility and Readonly -Visibility. The first 
one computes the visibility polygon in in-place manner with O(l) extra space. 
The second one assumes that the input array is read-only, and it uses 0(y/n) 
extra space. The time complexity of both the algorithms arc O(n). 

2.1 Inplace_Visibility 

We first describe an algorithm to report all the vertices of P that are visible from 
tt. Later, we show that the above algorithm can be easily modified to report the 
visibility polygon of tt. We start the algorithm by drawing a horizontal ray if 
through the point tt to its right that finds an edge eg = (pg,pg+\) of P intersected 
by 1$ first. Let q be the point of intersection pt and eg (see Figure 1(a)). 

We visit the vertices of P in counterclockwise order starting from pe+i- After 
visiting all the vertices of P, the visibility polygon of tt will be stored in consec- 
utive locations of the input array. Note that, at some point of time during the 
execution, a vertex may be visible to tt, but as the algorithm proceeds, it may 
not remain visible. 

Observation 1 Let pi and pj be a pair of visible vertices at an instant of time, 
where (i — 6) mod n < (j — 9) mod n. If both pi and pj become invisible during 
the further execution of the algorithm, then pj becomes invisible prior to pi . 

As in the classical algorithm for computing the visibility polygon [4], here also 
we store the vertices of P visible to tt in a stack. At the end of the execution, the 
content of the stack indicates the vertices of the visibility polygon of tt. Due to 



the constraint on space, here we maintain the stack in the array P itself. We use 
three index variables i, k and £, where k and I denote respectively the starting 
and ending indices of the stack at an instant of time, and i indicates the index of 
the current vertex of P under processing. We use two more workspaces # and \P 
that stores lqirP[£] and lqTrP[i] respectively, where t and i are as stated above. 
Note that the content of ^ can be from 0° to 360°. But & can contain negative 
angle. This happens when the traversal path along the boundary of the polygon 
crosses the ray Pi = nq an even number of times. 

The execution starts from the vertex P\6 + 1]. Initially we set k = I = + 1, 
i = 9+2, <P = & = 0. We process each vertex P[i], i = 6+2,9+3, . . . ,n, 1,2, . . . ,6 
in this order. While processing P[i], we compute \P = lqnP[i] as follows: 

• If P[i - 1] -> P[i] is an anticlockwise turn then \P = \P + LP[i - l]wP[i\; 

• Otherwise = V - LP\i - l]7rP[i]. 

Next, we compare $ with <P. Here one of the following three cases may arise. The 
appropriate actions in each case is explained below. At the end of the execution, 
the visibility polygon is available in P[k . . . , l\. 



Case 1: <P < 9 (see Figure 1(a)). 

Case 2: <P > & and the polygon makes a right turn at P[i] (see Figure 1(b)). 
Case 3: <P > & and the polygon makes a left turn at P[i] (see Figure 1(c)). 

In Case 1, P[i] is visible from the point n. We do the following: (i) push P[i] in 
the stack by setting I = [l + 1) mod n and placing P[i] in location P{(], (h) set 
<P = and (iii) increment i by setting i = (i + 1) mod n to process next vertex. 
In Case 2, we ignore P[i] and increment i (by setting i = (i + 1) mod n) until it 
encounters some vertex satisfying Case 1 (see Figure 2(a)). 
In Case 3, we pop elements from the stack by (i) setting <P = <P— LP\l\nP\l — 1] 
and (ii) decreasing £ by 1 at each step until one of the followings hold: 

Case 3.1: £ becomes less than k, i.e., the stack becomes empty. Here k 
and £ are reset to i (see Figure 2(b)). 

Case 3.2: P[£] and the current P[i] satisfy <P < & (see Figure 2(c)). Now, 
process P[i] as in Case 1. 

Case 3.3: P[£] and P[i] satisfy # > ^ and the line segments [P[i- 1], P[i}} 
and [P[£],P[£+ 1]] intersect in their interior (see Figure 2(d)). Here we 
need to proceed (by incrementing i by 1 at each step) until a vertex P[i'] is 
obtained that satisfies Case 1. Now we process P[i'] as in Case 1. 




Fig. 1: Processing of vertex P[i] 



Note that, in order to check Case 3.3, we need P[i — 1] and P[£ + 1] ; P[i — 1] may 
be lost during the execution of the algorithm. So we need to maintain P[i — 1] 
in a scalar location previous-vertex. P[£ + 1] is the last deleted vertex from the 
stack; it is available since no other vertex is inserted yet in the stack. 




(a)Caae2 (b) Case 3.1 («) Case 32 (d)Cnse3.3 



Fig. 2: Different cases that arise during the execution of our in-place algorithm 

Necessary modification required to obtain the entire visibility polygon 

In order to compute the entire visibility polygon, we need to modify the execution 
steps of Case 2 and Case 3 (where P[i] is not visible from ir) as follows: 

In Case 2, after computing the vertex i satisfying Case 1, we create a vertex 
4> at the point of intersection of the edge (pi-i,pi) and the line joining tt and 
P[£] (see Figure 2(a)). The point </> is pushed in the stack. Next, the vertex pi is 
processed as in Case 1. 

In Case 3, while backtracking (popping vertices from stack) by decrementing I 
by 1 at each step, we may arrive at one of the above three situations. 

In Case 3.1, q is a vertex of the visibility polygon. Thus, apart from resetting k 
and £, we need to insert q in the stack (see Figure 2(b)). Note that, during the 
entire execution, this case may appear at most once. 

In Case 3.2, a new vertex <j) is created at the point of intersection of the line 
joining [tt, P[i]) and the edge (P[£],P[£ + 1]) of P (see Figure 2(c)). Next, and 
P[i] are inserted in the stack. 

In Case 3.3, after computing P\i'\, a new vertex is created at the point of 
intersection of the line joining (n, P{£\) and the edge (P[i' — 1], P[i']) of P (see 
Figure 2(d)). Next, (j> and P[i'] are inserted in the stack. 

The pseudo-code of Algorithm 1 is given in the appendix. The correctness of the 
algorithm follows from Lemmas 3-7 of [5] and Lemmas 1 stated below. 

Lemma 1. The inplace maintenance of the stack in the same array P does not 
erase the polygonal vertices prior to its processing. 

Proof. Consider a push operation for maintaining the visibility polygonal vertices 
during the execution. Such a vertex may be either (i) a vertex of P or (ii) a point 
on an edge of P. When a polygonal vertex is pushed in the stack (Case (i)), i 
is immediately incremented. Thus, it does not erase any unprocessed polygonal 
vertex. In Case (ii), we need to mention that each pair of consecutive edges of 
the stack defines an edge of the visibility polygon [5]. While processing P[i], if a 
point 4* on the visibility polygonal edge (a, /3) is created, which is not a polygonal 
vertex, then both a and j3 are present in the stack. We pop a and push 4> in the 
stack. Thus the stack does not overlap unprocessed vertices. □ 



Apart from the array P, we need at most a constant number of scalar locations 
to run the algorithm. Its time complexity depends on the number of changes in 
the value of i and I. Since i is never decremented, i is modified n times. I is also 
modified 0(n) times since it is incremented at most n times, and the number 
of decrements of t is bounded by the number of its increments. If there is any 
reset in k during the execution, then £ is also changed. Thus changes in k does 
not affect the complexity of the algorithm. Thus, we have the following result: 

Theorem 1. The time complexity of the algorithm Inplace_Visibility is 0{n), 
and it uses 0(1) extra work-space. 

2.2 Readonly -Visibility 

Here, we assume that the vertices of the input polygon are given in a read-only 
array P in counterclockwise order, and show that the visibility polygon for a 
point 7r inside P can be computed in 0(n) time using 0(y/n) extra spaces. 

As in the earlier section, here also we draw a horizontal ray H from the point ir 
towards right that meets the edge eg = (pe,p$+i). We first partition the polygon 
into \y/n \ polygonal chains (polychains) ; each containing yjn vertices except the 
last one, which may contain fewer number of vertices. The first vertex of the first 
chain is pe+i- The j-th polychain will be referred to as Pj. We start processing 
from pe+i € Pi, and process the polychains in counterclockwise order. While 
processing Pj , the vertices of Pj are also processed in counterclockwise manner. 

Definition 1. A polychain Pj is called processed if all its vertices are processed. 
We maintain two arrays of integers: S of size \_y^n\ , and R of size 2 x \y/n \ . 

Let us assume that we have processed the chains Pi, ■ ■ ■ , Pj-i and next we want 
to process Pj. The array S contains the indices of the vertices in Pj. At an 
instant of time, if any part of the polychain Pj is visible from it, then R[l,j] 
and R[2,j] stores the indices of two vertices of P that blocks the visibility of 
Pj from its left and right sides at that instant of time. If the first (resp. last) 
vertex of Pj is visible to ir, then R[l, j] (resp. R[2,j]) stores that vertex itself. 
A zero entry in i?[l,j] and i?[2,j] indicate that Pj is entirely not visible to ir. 
Note that, if more than one part of Pj is visible from ir at an instant of time, 
then that information is not stored in the array R. From now onwards, by the 
term that a vertex is visible/invisible, we mean that it is visible/invisible to ir. 
The following two structural lemmas are crucial. 

Lemma 2. After the processing of Pj , let some vertices in the chain Pk (k < j) 
remains visible; the minimum and maximum indices of visible vertices in Pk be 
fk and l k , respectively. Now, if there exists any invisible vertex pp G P k with 
fk < (3 < lk, then the visibility of pp can only be obstructed by an edge of the 
chain Pk itself. 

Proof. If the visibility of pp is obstructed by some edge of an unvisited poly- 
chain, then it is not yet identified. So, we need to consider the case where pp is 
obstructed by some edge of a visited polychain P 7 , where l<7<fcor/c<7<j. 
In the former case, prior to obstructing pp, it must have obstructed p/ k . In the 
latter case, prior to obstructing pp, it must have obstructed pi k . In both the 
cases we have contradiction since both pf k and pi k are visible to ir. □ 



Lemma 3. During the processing of Pj, if it is observed that some visible ver- 
tices of Pk (k < j) becomes invisible, then all the vertices o/Pj_i, Pj-2, ■ ■ ■ , Pk+i 
(if any) becomes invisible. Moreover, if a vertex p a e Pj becomes invisible, then 
all the vertices having index greater than a in Pj are invisible. 

Our algorithm consists of two passes. In the first pass, we compute fj and lj 
for all the polychains Pj in anticlockwise order, and set R[l,j] and R[2,j] as 
described above. In the second pass, we consider each Pj and print the visible 
portions from it (if any). If R[l, j], R[2, j] ^ 0, then there exists at least one part 
of Pj that is visible from it. 

Pass 1: In this pass, we process Pj as in Subsection 2.1. We copy the indices of 
the vertices of Pj in array S, and process them in counterclockwise order. As in 
Section 2.1, during the execution, the visible portion of Pj is stored in the stack 
maintained at the beginning of S. An index variable I is used as the top pointer 
of the stack. We also use a variable \ that contains the index of the most recently 
visited polychain which is completely /partially visible at the current instant of 
time. At the beginning of processing Pj, \ is initialized with j — 1. 

While processing a vertex pe Pj, here also three cases may arise. The processing 
of different cases are same as in Subsection 2.1 except Case 2 and Case 3.1. 

In Case 2, if the last vertex p of Pj is not visible (i.e. lqirS{£] > Lq-Kp), then we 
store S[£] in a temporary variable a, We set R[2, j] = R{l,j + 1} = a, and copy 
Pj+i in S for the processing. 

In Case 3.1, the first task is to put the index of p in i?[l, j]. Next, we check 
whether the vertex p blocks the visibility of some already processed polychains by 
considering them in clockwise order starting from P x . By Lemma 3, a polychain 
Pk is considered if all the polychains {-P 7 ,fc < 7 < j} become invisible. While 
considering Pk, we may have the following three situations. Let the vertex p' be 
the clockwise neighbor of the vertex p e Pj, R[l, fc] — a and R[2, fc] = (3. 

Case 3.1.1 - The edge (p',p) does not intersect the lines [ir,p a ) and 
[it,pp): Here the entire visible portion of the polychain Pk remains visible, 
and we need not have to test the other visited polychains P 7 , 7 < fc. So, we 
push p on to the stack maintained in S, and consider the vertex next to p in 
the polygon P for processing. 

Case 3.1.2 - The edge (p',p) of P intersects the line [n:,pp) but does 
not intersect [ir,p a ): Here, we replace f3 by the index of p in R[2, fc]. Next, 
p is pushed in S and the vertex next to p is considered for processing. 

Case 3.1.3 - The edge (p' ,p) of P intersects both the lines [x,p a ) and 
[tt,pp): Here the entire polychain Pk becomes invisible. We set i?[l,fc] = 
R[2, fc] = 0, and consider the next visible polychain in clockwise order. The 
process continues until we arrive at either Case 3.1.1 or Case 3.1.2. 

At the end of processing Pj (if special Case 2, as mentioned before, doesn't arise), 
we put the index of the last vertex of Pj in R[2, j]. 

Pass 2: In this pass, we consider each polychain Pj for printing its visible 
portion from tt. We start from Pi, and proceed in counterclockwise order. If 



R[l, j], R[2, j] ^ for a polychain Pj, then it is fully/partially visible from tt. 
We compute the leftmost and rightmost points of Pj, say fj and lj that are 
visible from tt as follows. 

If R[l, j] points to the leftmost point of Pj, then fj = R[l,j]. Otherwise, we 
identify the edge e that is hit by the ray Wp^ first, where a'=R[l,j]. We 
compute <fr = point of intersection of e and TrpJ/ . 

Similarly, if R[2, j] points to the rightmost point of P i} then lj = R[2,j]. 
Otherwise, we identify the edge e that is hit by the ray Wp~/£ first, where 
/3' = R[2,j]. We compute ip = point of intersection of e and Trpf>. 

Next we copy the indices of all the vertices of Pj between <p and ip along with <p 
and ip in the array S in order. By Lemma 3, if there is any invisible vertex in 
S, it is blocked by some edge in Pj. So, we execute Inplace_Visibility on the 
array S. At the end of processing Pj, the content of the stack is printed. 
Theorem 2. The algorithm Readonly .Visibility correctly computes the vis- 
ibility polygon of P from the point tt. It needs 0{n) time and 0(y/n) work-space. 

Proof. The correctness of the algorithm Readonly .Visibility follows from the 
fact that for each polychain Pj, if its visibility is blocked from any/both side(s), 
then blocking vertices are correctly computed as mentioned in the correctness 
proof of the algorithm Inplace.Visibility and are stored in the array R. In 
Pass 2, the first and last visible vertices of Pj can be correctly computed from 
the content of R[l, j] and R[2, j]. The correctness of Pass 2 follows from Lemma 
2. We now analyze the complexity results of the algorithm. 

While executing Pass 1 on the polychain Pj , copying the polychain Pj in S needs 
0(y/n) time. During the processing of Pass 1, each vertex in Pj is inserted in and 
deleted from stack at most once. In the entire Pass 1, if a polychain becomes 
completely invisible, i.e., R[l, j], R[2, j] are set to 0, it never becomes visible. 
Note that, during the processing of each vertex in Pj, the visibility of at most 
one polychain is reduced by changing its j] field. Thus, the amortized time 
complexity for processing each vertex in Pass 1 is 0(1). While executing Pass 2 
for each Pi, we first spend 0(y / n) time for computing the first and last visible 
vertices <p and ip. The copying of the visible portion of Pi in S, and executing 
Inplace.Visibility on S needs another 0(y/n) time in the worst case. Since 
we have \y/n\ polychains, the result follows. □ 

3 Weak Visibility Polygon of an edge 

Given a polygon P and a line segment i = [p, q] in P, the weak visibility polygon 
of I, denoted by WVP(£), is a simple polygonal region R such that each point 
in R is visible from at least one point of I. As in the earlier section, here also 
we will assume that the vertices {pi,P2, ■ ■ ■ ,Pn} of the polygon P are given in a 
read-only array, called P, in anticlockwise order. We will propose an algorithm 
for computing the weak- visibility polygon of an edge e = (pi,Pi + i) of P, denoted 
by WVP(e) (see Figure 3). 

Before presenting the algorithm for computing WVP(e), let us consider first the 
following two simpler problems. 



LeftNon Visible (e, e'): Given two edges e and e' of a polygon P, compute the 

portion of e', from its left endpoint, which is not weakly-visible from e. 
RiGHTNONViSiBLE(e, e'): Given two edges e and e' of a polygon P, compute 
the portion of e', from its right endpoint, which is not weakly-visible from e. 
We will explain the method of solving the problem LEFTNONViSiBLE(e, e'). 
RightNonVisible(c, e') will be symmetric to that. Once we have the solution 
of these two problems, we will be able to compute the portion of e' which is 
weakly visible from e. In order to compute WVP(e), we need to compute the 
weak-visible portion of all the edges e' 7^ e of the polygon P. 

3.1 LEFTNONVlSIBLE(e, e') 

Let e = \pi,p i+ i] and e' = [pj,p j+ i], j > i. We first join [p i+1 ,p j+ i] (say I). We 
will use 9 and 9' to denote the end-points of L on e and e' respectively. At the 
end of the execution, it returns <P = left non- visible portion of e' from e. Initially 
9 = Pi+\ and 9' = Pj+\- The algorithm consists of three passes. 

Pass-1: In this pass, polychain 77 = \pj+i,pj + 2, . . . ,p n ,p\,P2, . . ■ ,Pi] is tra- 
versed in a counterclockwise manner starting from Pj+i- If a vertex p £ 77 is 
observed which is to the right of L in its present position, then 9' is moved to 
the point of intersection of the line containing e' and the line (pi+i,p); 9 remains 
fixed at Pi+\. If 9' is outside e', then e' is not visible from e, and the procedure 
returns <P = e' . Otherwise, L is defined by p i+ i and a vertex p T € II (see Figure 
4(a)). A temporary variable r is used to remember p T . 

Pass-2: In this pass, we simultaneously traverse 77 from pi in clockwise order up 
to the vertex p T , and 77' from Pi+i to pj in counterclockwise order. The method 
of traversal is explained in Process-1 and Process-2, stated below. We use three 
index variables fc, £ and m. Initially, we set k = i+l, £ = i and m = r. 

Process-1: We traverse 77 in anticlockwise direction using the index variable k. 
At each move one of the following events may be observed: 

• k = j. Pass 2 stops, and it returns <S> = [pj+i, 9']. 

• k 7^ j and the edge (pk-i,Pk) does not intersect L. Here k is incremented 
(modn) to process the next vertex of 77. 

• k 7^ j and the edge (pk-i,Pk) intersects L above p m . Here, e' is not weakly 
visible to e (see Figure 4(a)), The procedure returns <P — e! . 

• k 7^ j and the edge (pk-i,Pk) intersects I below p m . Here, we update L with 
the line joining (p m ,Pk), and update 9 (resp. 9') by the point of intersection 
of L and e (resp. e') (see Figure 4(b)). Next we execute Process-2. 



Fig. 3: Demonstration of (a) Weak visibility polygon, (b) minimum link path 




Fig. 4: Algorithm LeftNonVisible 



Process-2: We traverse 77' in clockwise direction using the index variable £. At 
each move the following events may happen: 

• the edge (pe+i,pg) does not intersect L. Here, I is decremented (modn) to 
process the next vertex of 77'. 

• the edge (pe+\,pi) intersects L below pk- Here, e' is not weakly visible to e 
(see Figure 4(c)). The procedure returns <P = e'. 

• the edge (pt+i,pe) intersects L above pk- Here, we update the L by the line 
joining pt and p^. We also set 6' = ir and m = £, and then switch to executing 
the Process- 1. 

Pass-3: Note that after executing Process 2, L may intersect 77 below pk (see 
Figure 4(e)). Thus the weak-visibility of e' from e may be lost. So, after the 
successful finish of Pass-2, we execute Pass-3 to check for possible intersection of 
L and the edges of 77'. If such an intersection is observed, e' is not weakly visible 
from e; otherwise, the final position of 6' determines the non-visible portion 
<P = \pj+i,0'] of e' from e (see Figure 4(f)). 

Lemma 4. Algorithm LeftNonVisible correctly computes the invisible por- 
tion of e' from its left end-point, and it needs 0(n) time and 0(1) space. 

Proof. Initially, we draw the line segment L = [pi + i,pj + {\. Now, three cases may 
arise: (i) the entire segment L lies inside the polygon, (ii) the entire segment L 
lies outside the polygon, and (iii) L spans both inside and outside the polygon. 
In case (i), e' starting from Pj+i is visible. In other two cases, there may exist 
situations where e' is/is not weakly visible from e. We now show that in both 
cases our algorithm correctly computes the left non- visible portion of e' . 



In Pass-1, we compute the invisible portion of e' from left by traversing the left 
chain 77 = [pj+\ ~ Pi] (assuming the right chain is absent). At each step of 
traversal, if the current edge obstructs L, then L is shifted towards right along 
e', and in the further scanning of the chain L is never shifted towards left. 

Now, we need to consider the right chain 77'. The visibility may be blocked by 
some vertex of 77' at the current position of L. So, we traverse the right chain 
from p i+ i to pj. As soon as a blocking is observed, L is modified as mentioned 
in the algorithm. However, in the new position of 7, a portion of it becomes 
closer to the left chain 77; thus, it may again be obstructed by 77. So, we need 
to traverse 77. This alternating process may continue until pj is reached along 
the right chain. It needs to be mentioned that during this traversal, we have 
not noticed whether L is obstructed by some edge which is already visited. So, 
finally we execute Pass-3 to check this. 

The time complexity follows from the fact that 77 is traversed three times and 
77' is traversed once only. During the traversal of 77 (resp. 77') its each vertex 
is visited at most once. □ 

Theorem 3. The weak visibility polygon of an edge of P can be computed in 
0(n 2 ) time using 0(1) extra space, where the vertices of P are given in a read- 
only array. 

4 Minimum link path between a pair of points 

Given a polygon P and a pair of points s and t, the minimum link path between s 
and t, denoted by MLP(s, t), is a polygonal chain from s to t where the number 
of edges in the chain is minimum among all other polygonal chains connecting s 
and t (Figure 3(b)). We propose an algorithm for computing MLP(s, t) assuming 
that the vertices of P are given in a readonly array in anticlockwise order. 

A classic way to compute MLP(s, t) is as follows: (i) Compute the visibility 
polygon of s, called Q. If t € Q, then s and t are straight-line visible. Otherwise 
(ii) identify the edge \ 01 Q such that the sub-polygons of P lying on one side of 
X contains s and that on the other side contains t. Now (iii) compute the weak 
visibility polygon of \ m the sub-polygon containing t. If again it contains t, 
then the process stops; otherwise iterate steps (ii) and (iii). The time complexity 
is 0(kn 2 ) and it needs 0(1) extra space (see Theorem 3), where k is the number 
of segments in MLP(s, t). 

We now describe an algorithm for this problem that runs in 0(kn) time and 
0(1) work-space. It also executes k iterations to report the k segments of the 
path. We use a variable it to store the point such that the path from s to it is 
already reported. Initially tt stores s. In each iteration we execute the procedures 
LeftNon Visible (e, e') described in Section 4. It is tailored such that it can 
work even if e and/or e' is/are point (s) inside P, and returns two parameters <P 
and L, where is the left non-visible portion of e from e' = t and 7 is a line 
such that nothing to the left of 7 is visible to e' = t. If e is a point, then we draw 
a horizontal line segment [a, b] inside P that passes through e and the points 
a, b lying on the boundary of P. Similarly, the line segment [c, d] is defined for 



(a) (b) (c) (d) (e) 

Fig. 5: Computation of minimum link path 

e' if it is a point. Initially, e and e' are the points s and t, respectively. Here 
77 is defined as an anticlockwise polychain from a to c and 77' is a clockwise 
polychain from b to d. From next iteration onwards e = [a, b] will be an edge of 
the corresponding polygon, but e' remains equal to t in all the iterations. 

At each step, if <P ^ e and one of the end-points of L is equal to t, then t is 
visible from some point of e. We choose the left-most point q of and report 
the last two edges [n,q] and [q,t] of the Minimum_Link_Path. 
If $ = e, then no part of e is visible from t; We execute the procedure Compute- 
First-Link(<£, L), stated below to compute an edge \P of the weak- visibility poly- 
gon of e such that s and t are in two different sides of W. 

Observation 2 (i) Every point of & is weakly visible from e and no point in 
the proper interior of the sub-polygon of P containing t is visible from e. 
(ii) Extension of & intersects e (at a point, say q). 

We report the link [tt, q]; reset tt — q, and execute the next iteration with e = \P. 
4.1 Compute-First-Link 

Here we need to handle the following two cases: (i) L does not intersect the edge 
e but contains the point t, and (ii) L does not contain the point t. 

Case (i) Let L be defined by a vertex pi e 77 and t (see Figure 5). We com- 
pute the right non- visible portion of e = [a, b] from the point pi. We redefine 
L = [b,pi] and alternately visit the vertices of the polychain 77' from d to b in 
anticlockwise order and the polychain 77 from pi to a in clockwise order using 
the procedures Process-1 and Process-2 stated below. Note that, [c, d] is a 
chord of P containing t, as stated earlier. 

Process-1: 

• If no edge of 77' intersects L above pi then it returns W = [pi,p], where p is 
the point of intersection of L and 77' (see Figure 5(a)). 

• If an edge (pk,Pk+i) of 77' intersects 7, then L is modified by [pi,pk), and 
the end-point of L on e moves towards a. For each encounter of an edge of 
77' intersecting L, L is modified accordingly. The traversal continues along 
77' until b is reached (Figure 5(b)) or an edge (pe,pe+i) € 77' is encountered 
such that the redefined L docs not intersect e (Figure 5(c)). Now Process-2 
is invoked. 



Process-2: 

We redefine L = [p e , a] and start traversing 77 from pi towards a to compute 
the left non-visible portion of e from the point p(. As in Process-1, for each 
encounter of an edge (p m ,p m+1 ) e 77 intersecting L, the end-point of L on e 
moves towards b. Finally 

• If a is reached and L does not leave e, then report = [p m +i, p], where p is 
the point of intersection of L and the polychain 77' (Figure 5(c)). 

• If the end-point of L on e goes beyond b for an edge (p m ,p m +i) € 77, then 
L is redefined as [p m+ i,b] and Process-1 is invoked. 

Case (ii) Let L be defined by two vertices pi e 77 and pj e 77'. As in the 
algorithm LeftNonVisible here also we traverse from c to pi along 77 to get 
a vertex that intersects [pj,Pi]. If no such vertex is found and L intersects 77 
at a point p, then = [p,Pj] is returned (see Figure 5(d)). If such a vertex pi> 
is found, then L is redefined by pj and pi< , and the traversal from pj starts to 
get a vertex in 77' that intersects [pj,pi'} (see Figure 5(e)). This type alternate 
traversal in 77 and 77' ultimately defines a line L by two vertices p a € 77 and 
pp £ 77' such that there exists no edge (pk,Pk+i) in 77 with k > i that intersects 
L. Now if L does not intersect e, the situation is similar to Case (i). Otherwise, 
the traversal along 77' continues from pp until it is obstructed by a vertex in 77. 
Then the traversal starts from p a e 77. This type of alternate traversal in 77 and 
77' continues as in Case (i) to have a chord of P defined by its two vertices. 

Theorem 4. The proposed algorithm correctly computes the minimum link path 
between s and t in time 0(kn), where P is given in a read-only array, and k is 
the size of the output. 

Proof. The above scheme reports an edge i] of the visibility polygon of e such 
that s and t are in different sides of it. The minimality of the path length follows 
from the fact that -q satisfies Observation 2. The above scheme needs 0(n) time 
since it needs visiting all the vertices of 77 and 77' at most a constant number 
of times. Since the length of the minimum link path is k, the result follows. 
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Appendix 



Algorithm 1: InplaceJVisibility 



Input: A simple polygon in the array P[] and a point 7r inside the polygon 
Output: Visibility polygon of 7r 

1 Draw a horizontal line H from tt towards right and find the first intersection point q with 
the polygon P. Let q £ the edge (p# , Pe+i )■; 

2 Initialize: k — I — Q + 1 •. i — 9 + 2; previous _vertex — P[9 + 1]; 

3 while i ^ 8 do 

4 (* Case 1: *); 

5 if LqnP\l\ < LqitP\i\ then 

6 i Push_and_Proceed(i, £); 

7 Case 2: *); 

8 else if ^qnP[£] > AqtrP[i] and the polygon makes a right turn at P[i] then 
while (lqnP[e] > /qTrP[i])/\(i # 9) do 

j Only_Procced(i) 
if Lq-nP\t\ < LqnP\i\ then 

<f> — intersection point of edge (pi-i,pt) and the half-line [7r,P[£]); 
Push_Phi(0, 1); 
PUSH_and_PROCEED(i, £); 
(* Case 3: *); 
else 

while !(Cose3.1 V Case3.2 V Case3.3) do 

I i = i-l; 
(* Case 3.1: *); 
if £ < k then 
i k = i - 1; P[k] = q; e = i\ 
I Only_Procced(i) 
(* Case 3.2: *); 

else if lqnP[e] < /qirP[i] then 

4> — intersection point of edge (P[^],P[j 
Push_Phi(0,£); 
PUSH_and_PROCEED(i, I); 
(* Case 3.3: (/<jttP[£] > lqirP[i\) A ([P[i - 
intersect) *); 
else 

while !{{LqnP[t] < Z<?7rP[i]) = 0)) do 

i Only_Procccd(i) 
if i ^ 6 then 

tp — intersection point of the edge (previous-vertex, P[i]) and the half-line 

, i Push_Phi(0,£); 

PUSH_and_PROCEED(i, £) 
of the array; 



I + 1]) and the half-line [tt, P[i]); 



l],P[i]] and [P[£],P[^ + 1]] properly 



35 

36 Report P[k . 

37 end. 

38 — 



39 Procedure PUSH_and_PROCEED(i : I); 

40 £ — (£ + l)mod n; previous .vertex — P[i\; 

41 swap(P[z], P[£]); i = (i + l)mod n: 

42 end. 

43 

44 Procedure Only_Proceed(i) ; 

45 previous .vertex — P[i\; i — (i + l)mod n; 

46 end. 

47 

48 Procedure Push_Phi(0, €)\ 

49 £ = (£ + l)mod n; P[£] = <p; 

50 end. 



Algorithm 2: Modificd-Case-3.1 



Workspace: An array Ft of size 2 X [y^n] ■> an< ^ array S of size LV^-J 
R[l,j] and R[2,j] contain the blocking vertices of Pj from its two sides; 

While processing Pj , it is copied in S and processed there. The array P remains unaltered; 
Variables used; 

I: top pointer of stack maintained at the beginning of S; 

X- index of the poly-chain of maximum index that is totally /partially visible: 

i: index of the current vertex in S; 

(* Modified Case 3.1: Here I is reached to *); 

while x] - A R[2, x] — 0) V ((previous -vertex, S[i]) intersects both the line [ir,p a ) 

and [7r : pj3)) V (x 7^ 0) do 

if (previous ^vertex, S[i]) intersects both the line [ir,p a ) and [ir,pp) then 
| (R[hx\ = 0; R[2,x] =0; 

x = x - i; 

if x ¥= then 

if (prev ious-vertex, S[i]) intersects none of the line [-7T,p a ) and [tr^pp) then 

i previous .vertex — S[i}; I = t + 1; S[£] — S[i]; i = i + 1; 
else 

j H[2,x] = »; 

! previous ^vertex — S[i]\ i = I + 1; 5[^] — z = z + 1; 



Algorithm 3: Pass 2 of Algorithm Readonly -Visibility 



j = i; 

while 
if 



16 
17 



j ^ No_of -Partition do 
R[l, j] ^ A i?[2, j] then 
a = K[l,j]; /3 = R[2,j]; 
if P[ct] g Pj then 
1 Find the edge e of Pj which is intersected first by the line [-7T, P[a); 
! <f> — intersection of e and [-7T, P[a]) 
else 

I 4> = P[<*} 

if P[0] g Pj then 

j Find the edge e of Pj which is intersected first by the line [tt, P[/3); 

| ip — intersection of e and [tt, P[{3]) 
else 

I ■<!> = P[P] 

Copy to i/? of the polychain S', 

Execute Algorithm INPLACE.VISIBILITY on S to report the visible portion 

= j + i 



Algorithm 4: LeftNonVisible(^^ +1 , VjVj+v 



Input: Two edges pipi+i and PjPj+i of a simple polygon P 

Output: The left portion of pjPj+i which is not weakly- visible from PiPi+i 

PASS1Q; 

r = PASS2Q: 

if r — — 1 then 

Report nothing is visible; 

else 

| PASSSQ; 
end. 



Procedure PASSlQ; 

L = (Pj+i,Pi+i)\ 
t = {j + 2)mod n; 
while £ 7^ i do 

if p± is to the right of L then 
i UpdateLine(p £ , Pi+i ); 
r = * ; 
t = (t + l)mod n; 

end. 



19 Procedure PASS2(); 

20 k — i + 1; I — i; m = r; 

21 r — 1; 

22 while r/OAr/- 1 do 

23 if r — 1 then 

24 l Processl(); 

25 else if r — 2 then 

26 j Process2(): 

27 Return r; 

28 end. 



30 Procedure ProcesslQ 

31 if fc — j then 

32 | r — 0; 

33 else 

34 if (pfc_iPfe) does not intersect line L then 

35 j k — (k + l)mod n; 

36 j r — 1; 

37 else if (pfc_ipfc) intersects below Pm then 

38 | UPDATELlNEfp™., p k ); 

39 j r — 2; 

40 else 

41 | r = — 1; 

42 end. 

43 

44 Procedure Process2() 

45 if I > r then 

46 if p^+ip^ does noi intersect L then 

47 I £=(£-1) mod n; 

48 ! r — 2; 
else if pi+\pi intersects L below p k then 

r = -1; 

else 

UpdateLine(pz , pfc); 
m = £ : 
r = 1 ; 



55 else 

56 | r — 0; 

57 end. 



59 Procedure PASS3{) 

60 t- (j -f 2)mod n; 

61 while (t / z) /\(r / -1) do 

62 if p± is to the right of L then 

63 | r = — 1; 

64 t = (t + l)mod n; 

65 £ — (i + 2)mod n; 

66 while (£ / j) /\(r / -1) do 

67 if p t is to the left of L then 

68 j r = — 1; 

69 t = (t + n; 

70 if r — —1 then 

71 | Report nothing is visible; 

72 else 

73 | Report (pj is not visible; 

74 end. 



